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We propose a scheme to eliminate the effect of non-nearest-neighbor qubits in preparing cluster 
state with double-dot molecules. As the interaction Hamiltonians between qubits are Ising-model 
and mutually commute, we can get positive and negative effective interactions between qubits to 
cancel the effect of non-nearest-neighbor qubits by properly changing the electron charge states of 
each quantum dot molecule. The total time spent on the present multi-step cluster state prepara- 
tion scheme is only doubled for one-dimensional qubit chain and tripled for two-dimensional qubit 
array comparing with the time spent on previous protocol leaving out the non-nearest-neighbor 
^— ^ ' interactions. 
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Semiconductor quantum dot (QD) is believed to be one of the most promising systems for quantum information 
processing [1, 2]. Recently, quantum molecule, formed by coupled double quantum dots, has attracted many interests. 
A lot of experimental [3, 4, 5, 6] and theoretical [7, 8, 9] work has be done on the two-electron states of the double-dot 
. molecule. It is argued that by encoding in singlet and triplet states of quantum molecule, qubits can be protected 
from low-frequency noise and the dominant source of decoherence from hyperfinc interactions can be suppressed 
Q\ , [3, 10, 11, 12, 13, 14]. Universal quantum gates and a fault-tolerant architecture for quantum computation have been 
■ proposed for these qubits encoding in the two-electron states of quantum molecule [8, 15]. In the reference [16], we 
proposed a one-step scheme to prepare large cluster states with QD molecules, which can act as a general source for 
one-way quantum computation [17]. As in most of the solid-state quantum computation protocols, only interactions 
between neighboring qubits are considered in this scheme. However, all these protocols directly explore the long range 
Coulomb interaction of the electrons between different molecules [8, 15, 16]. There are unwanted small interactions 
between non-nearest-neighbor molecules. Generally, this kind of unwanted interactions between non-nearest-neighbor 
qubits will act as an additional decoherence source and affect the proposed quantum processing fidelity, for example, 
the fidelity of prepared states or quantum gates. It is well known that there is a decoupling and refocusing technique in 
standard nuclear-magnetic-resonance field [18, 19]. Many works explores this idea to simulate quantum dynamics [20], 
construct universal quantum computation with some Hamiltonian forms [21]. Here we focused on a particular physical 
system of double-dot quantum molecules, and propose an efficient scheme to eliminate the effect of the interactions 
between non-nearest-neighbor qubits in the preparation of cluster state with QD molecules. By properly changing the 
electron charge states in each QD molecule, the effective interactions between qubits can be switched as positive, zero 
and negative. As the interaction Hamiltonians between qubits are all Ising-model and mutually commute, the effect 
of non-nearest-neighbor qubits interactions can be efficiently canceled through a multi-step preparation process. In 
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addition, the total time for the present multi-step cluster state preparation scheme is only several times comparing 
■ - - i with the time of previous protocol neglecting the non- nearest-neighbor interactions. 

For each double-dot quantum molecule, the notation (n u , rii) can be used to indicate the number of electrons in the 
upper and lower QDs. Define a bias parameter A to represent the potential offset between the three possible charge 
states of each molecule (0,2), (1, 1) and (2,0). We can adiabatically sweep A of each molecule by tuning gate-bias 
voltages of each molecule or applying an external electrical field. The charge states (0, 2), (1, 1) and (2, 0) respectively 
correspond to the case of A = E c ,0 and — E c as shown in Fig.l [3, 7]. 

For the charge state (1, 1), there are four spin states: \S) = (|U) - |lt))/V2, \T) = (|||) + |1T»/V2, \T+) = ITT), 
and |T_) = \H). The two triplet states \T±) can be largely separated from the singlet state \S) and triplet state 
\T) by Zeeman splitting, which enables us to exclude them in the following [8]. The qubit can be encoded in the 
near-degenerate (1, 1) charge states \T) = |0) and \S) = which are separated due to weak tunneling between the 
two dots of each molecule. As a result of Pauli blockade, only the singlet state \S) can be tuned into charge state 
(0, 2) or (2, 0) by sweeping A. The triplet state |T) will be kept in the charge state (1, 1) in the processing of sweeping 
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A. 

As the Coulomb interaction between the two electrons inside each molecule can be excluded, we only consider the 
interaction between electrons of different molecules [15, 22]. For any two quantum molecules, different charge states 
can result in different Coulomb interactions as shown in Fig. 2. When the bias voltage of two molecules are both set on 
A = 0, these two molecules can be both in the charge state (1, 1) (see Fig. 2(a)). In this case, the Coulomb interaction 
energy between them can be directly written as: 

E (g,b) = — (^-+ ^ V (1) 
V ' 47re V b Va 2 + b 2 J 1 ' 

Here and in the following text, e is dielectric constant of GaAs, a is the distance between the two QDs in one 
molecule and b is the distance between the two molecules. When one of two molecules remains in the charge state 
(1, 1) by keeping its A = and the other molecule is changed to other charge state (0,2) or (2,0) by sweeping its 
A to E c or — E c , the interaction between them can remain in the form Eo(a,b) as shown in Fig. 2(b). When the 
two molecules are both in the charge state (0,2) or (2,0) by tuning their bias parameter A, the interaction energy 

between them can then be written as E\{a,b) = j^^f- as shown in Fig. 2(c). As only the singlet state \S) can be 
tuned into charge state (0, 2) or (2, 0) and the triplet state \T) will be kept in the charge state (1, 1) by sweeping A, 
this interaction Hamiltonian between two neighboring qubits can be written in the basis \TT), \TS), \ST) and \SS) 
as [16]: 

H+(a, b) = diag{£ (a, b),E (a, 6), E {a, b), E 1 {a, b)} . (2) 

Taking the interaction energy E (a,b) as zero point of the effective interaction energy between two qubits, we can 
write the interaction Hamiltonian as: 

H+(a,b) = H (a,b) + H' + (a,b), (3) 

here we define: 

H (a,b) = E (a,b) x diag{l, 1, 1, 1} ; (4) 



H' + (a,b) = diag{0,0,0,£+(a,6)}, (5) 

where E + (a,b) = Ei(a,b) — E a (a,b) = ^+v 2 ) * s P osrt i vc - When one of two molecules is in the charge 

state (0,2) and the other one is in the charge state (2,0) as shown in Fig. 2(d), the interaction energy between them 
can be written as E 2 (a,b) = j^-^===. Similarly, we can also get effective interaction in the basis \TT), \TS), \ST) 
and \SS) as: 

H-(a,b) = dmg{E (a,b),E {a,b),E (a,b),E 2 (a,b)} 

= H (a,b) + H'_(a,b), (6) 

where 

H' + (a,b) = diag {0,0,0, (a, 6)}, (7) 

and E-(a,b) = E^ia^b) — Eo(a,b) = —E + (a 1 b) is negative. 

The effective interaction between two non-neighboring qubits can also be written in the above forms, except that 
the parameter b is replaced by the distance between the two non-neighboring molecules. Therefore, we can set the 
effective interaction between any two quantum molecules as positive, zero and negative by properly changing the 
charge state of each molecule. The property enables us to cancel the effect of non-nearest-neighbor qubits interactions 
with a multi-step preparation process for cluster state of QD molecules. 

Firstly, all qubits are initialized in the states (\S) + \T))/\/2. In this case, all molecules are in the (1, 1) charge state 
and the effective interaction between any two qubits is closed. Sweeping A of all molecules by tuning gate-bias voltages 
of each molecule or applying a globe external electrical field, the qubits in the IS*) state can be changed into (0,2) 
charge state. With this operation, Ising-model interactions as Eq.2 are switched on between all quantum molecules. 
As in our previous scheme [16], if there are only interactions between nearest-neighbor quantum molecules, we can 
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prepare cluster state simply with this operation. However, there are small interactions between non-nearest-neighbor 
QD molecules. For example, there is interaction between qubit i and qubit i + k in the following form: 

H + (a,kb) = E (a, \k\b) x diag{l, 1,1,1} + 
diag{0,0,0,£+(a,|fc|&)} 
= H (a,kb) + H' + (a,kb). (8) 

Here kb represents the distance between molecule i and i + fc; k = ±1, ±2, • • • ; when k < 0(k > 0), i + k correspond 
to the qubits of left (right) side neighbor. We describe the effect of non-nearest- neighbor qubits by the ratio of all the 
non-nearest-neighbor interactions to the nearest-neighbor interaction (based on the following discussion, the effect of 
Ho(a, kb) can be ignored): 

R{a,b) = Y,E + {a,kb)/E+{a,b)- (9) 

fe=2 

For a typical value of a and b = 10a, the ratio R(a, 10a) s=s 20%. Therefore, we need to include interactions between 
non-nearest-neighbor qubits to generate a cluster state with higher fidelity. 

Including the non-nearest-neighbor interactions, the total interaction Hamiltonian can be written as: 

+ 0O / +0O \ 

H= Y, [H? +1 (a,b) + YH? +k (a,kb)\. (10) 

i=-oo V k=2 J 

The superscripts in Eq.10 indicate the two qubits which the interactions acting on. It is noted that the interaction 
Hamiltonians between any two molecules are Ising-model and mutually commute, so the order of application of the 
time evolution operations does not matter. We can then describe the system time evolution as follows: 



exp 



I E [Hr +1 (a,b) + J2Hr +k (a,kb) 




Here 

U hl+k (t) = exp 



itH + (a, kb) 

h 



= e^)diag{l,l,l,ex P ( ig+ y^ )}, (12) 



where 



e*<"> = exp (^«) , (13) 

is a constant phase, which only affect the trivial global phase of final state. That means we can ignore the Ho(a, kb) 
part in H±(a, kb). For briefness, we use H±(a, kb) = H±(a, kb) 

If there are only interactions between nearest-neighbor quantum molecules, the evolution operator includes simply 
the term n£^oo U l ' t+1 , with which the system can evolve into cluster state for a proper time t = t [16]. In order to 
eliminate the effect of non-nearest-neighbor qubit interaction, we need to cancel the operation nS^oo k=2 U 1,l+k (The 
constant phases caused by each Ho(a, kb) can be added together as a global phase, which has no effect to the final 
state fidelity). As shown above, we can get negative Ising-model interaction H_(a, kb) = — H + (a, kb) by changing the 
molecule charge state. If the system successively evolves under H + (a, kb) and H_(a, kb) for a same time interval, the 
system will return to its origin states. Thus we can split the cluster state preparation into several steps to eliminate 
the effect of non- nearest-neighbor qubits: after the first step operation 17(f), we can properly switch on negative 
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interaction (a, kb) to cancel the effect of the operator Jlit^oo k=2 U ht+k . The number of the steps depends on how 
many non-nearest-ncighbor qubit interactions we want to cancel. 

For example, to eliminate the effects of the next-nearest-neighbor qubits, we can use a three-step process as given 
in Fig. 3 and Tabled to generate cluster state. Firstly, after initialization, wc tune bias parameter A of each molecule 
to move qubits i + 4n, i + An + 1 (the number n = 0, ±1, ±2, • • • ) into the charge state (0, 2) and qubits i + An + 2, 
i + An + 3 into the charge state (2,0) for a time interval of ti = to/2 (Fig. 3(a)). Here, to is the time we need to 
generate cluster state without considering the interactions between non-nearest-neighbor qubits. Then we move the 
qubits i + An, i + An + 3 into the charge state (0, 2) and qubits i + An + 1, i + An + 2 into (2, 0) for a time interval 
of t2 = to/2 (Fig. 3(b)). Finally, we can move all molecules into the charge state (0,2) for a time interval of t% = to 
(Fig. 3(c)). Through these three steps causing a time of 2t , we can prepare a cluster state where the effect of all 
the (An + 2)-nearest-neighbor qubits (n = 0, ±1,±2, ■■■ ), including next-nearest-neighbor qubits (n = 0,-1), are 
eliminated. Here (An + 2)-nearest-neighbor means the qubit separated An + 2 sites from qubits i. By indicating 
the effective interaction energies of each step in Tabled, we can understand this preparation process more clearly. 
Although the interactions between qubit i and qubit i + An are doubled, the ratio of the residual non-nearest-neighbor 
interaction to the nearest-neighbor interaction is reduced to the value R' through this three-step process: 

Rl = En=i(E+(a, (An ± 1)6) + 2E+(a,Anb) 
E+(a,b) 

w 9%, (b = 10a). (14) 

To eliminate more non-nearest-neighbor interactions, we can use a more-step process as shown in Table. II. To cancel 
the effect of (mn + 2)-nearest-neighbor to (mn + m — 2)-nearest-neighbor qubits (to is an integer lager than 3, and 
n = 0, ±1, ±2, • • • ), firstly, we move the qubits i + mn, i + mn + 1 into (0, 2) state and the others into (2, 0) state 
for a time interval of t\ = to /A (to is the same as before). Secondly, we move the qubits i + mn — 1, i + mn into 
(0, 2) state and the others into (2, 0) state for a time interval of ti = to/ A. Thirdly, we move the qubits i + mn — 2, 
i + mn — 1 into (0, 2) state and the others into (2, 0) state for a time interval of t 3 = to /A. Similarly shift the charge 
state sequence one qubit to the left for a time interval of to /A as above to times. After that, we move all the molecules 
into the charge state (0,2) for a time interval of t m +i = (8 — m)to/A. All the effective interactions between qubit 
i and qubit i + k in this process are shown in Table. II. It is demonstrated that after the above to + 1 steps, only 
the interactions between qubit i and i + nm ± 1, i + nm are left (the interactions between qubit i and i + nm are 
doubled). Thus we can prepare a cluster state through this process, where the effect of (mn + 2) -nearest-neighbor to 
(mn + m — 2)-nearest-neighbor qubits is eliminated. To satisfy t. m+ i = (8 — m)to/A > (the time interval of the last 
step should be positive), to must be no larger than 8. The ratio of the residual non-nearest- neighbor interactions to 
the nearest-neighbor interaction is reduced to the value R": 

El=i(E+(a, (mn ± 1)6) + 2E + (a, mnb)) 
E+(a,b) 

w 1%, (b= 10a, to = 8). (15) 

The total time we need to prepare cluster state in this multi-step process is only doubled. It is noted that this above 
process only can be applicable for the case 4 < m < 8. When m = 3, it corresponds to the case of eliminating the 
effect of next-nearest-neighbor qubits. There are still some non-nearest-neighbor interactions left in this cluster state 
generation process. However, as we can turn on or off the interaction H + (a, kb) or H_(a, kb) between any two qubits, 
the remaining non-nearest-interactions can also be eliminate by a more-step process, where the negative interactions 
H-(a,kb) are properly switched on and off to cancel the effect of non- nearest-neighbor interactions H + (a,kb). It 
is noted that the present idea can only be applicable when the interaction Hamiltonians between qubits mutually 
commute. 

In the above section, only one-dimensional qubit chain is considered. We can expand the above idea to a two- 
dimensional qubit array, where the two dots in each molecule can be fabricated in different layers of a bilayer system 
as shown in Fig. 4 [23, 24, 25, 26]. Each molecule can be transferred among (0,2), (1, 1) and (2,0) charge states by 
tuning the voltage of top gates and back gates. For example, we can eliminate the effect of the next-nearest-neighbor 
(diagonal) qubits interactions in a three-step cluster state preparation process. Firstly, we move qubits (i + 2m, j + 2n) 
to (0, 2) state, qubits (i + 2m + 1, j + 2n+ 1) to (2, 0) state and the other qubits to (1, 1) state for a time interval of t 
(Here to, n — 0, ±1, ±2, • • • ). Secondly, we move qubits (i + 2m, j + 2n + 1) to (0, 2) state, qubits (i + 2m + 1, j + 2n) 
to (2,0) state and the others to (1, 1) state for a time interval of to- Finally, we add a negative voltage on the top 
gates to move all the initialized qubits to (0, 2) charge state for a time interval of to- With these three steps causing 
a total time of 2>to, we can prepare a two-dimensional cluster state while the next- nearest-neighbor (diagonal) qubits 
interactions are eliminated. Similarly, with a more-step process, other non-nearest-neighbor interactions can also be 
eliminated by properly switching on the positive or negative interactions between each two qubits. 
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In real experiment, there are unavoidable randomness in the the distances between quantum dots (the value of a 
and b). This randomness will greatly affect the fidelity of final state, especially when the interactions between non- 
neighboring qubits are included. As in the present protocol, the interactions of non-neighboring qubits are canceled 
by switched to negative and positive value. The randomness in the distances between non-neighboring qubits is also 
eliminated. 

For a wide class of so-called "toy models" in quantum computation, most researches focus on the idea rather than its 
quantitative development in real systems. Nevertheless, the importance and efficiency of those simplified models are 
still very high, since they usually provide physical ground for more involved theories. Moreover, the lack of pronounced 
results in up-to-date experimental efforts in realization of even prototypical quantum registers, composed of tens of 
qubits, forces us to consider and to analyze all proposed frameworks of quantum computations. Our results may 
stimulate further theoretical investigations of perspectives of one-way quantum computations for which the existence 
of reliable mechanism of cluster state preparation is crucial. And it be interesting and useful for researchers dealing 
with this area of investigations. 

In conclusion, we propose a multi-step cluster state generation protocol for double-dot molecules to eliminate the 
effect of the interactions between non- nearest- neighbor qubits. As the interaction Hamiltonians between qubits are 
Ising-model and mutually commute, the order of application of the time evolution operations does not matter. We 
can get positive and negative effective interactions between qubits to cancel the effect of non-nearest-neighbor qubits 
by properly tuning the electron charge states of each QD molecule. The total time for the present multi-step cluster 
state preparation scheme is only doubled for one-dimensional qubit chain and tripled for two-dimensional qubit array 
compared with the time of previous protocol neglecting the non-nearest-neighbor interactions. 

This work was funded by National Fundamental Research Program, the innovation funds from Chinese Academy 
of Sciences and National Natural Science Foundation of China (Grant No. 60121503 and No. 10604052). 
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TABLE I: Interaction energies between qubit j and qubit j+k in the process of eliminating the next-nearest-neighbor interactions 
in preparation of cluster state(qubit j is any qubit in the molecule chain). The rows begin with (a), (b) and (c) are corresponding 
to the Fig.3(a), Fig. 3(b) and Fig.3(c) in Fig. 3, respectively. And the numbers in last four lines indicate the effective intense of 
the interaction between qubit j and qubit j + fc(Units: E+(a, kb)). 
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TABLE II: Interaction energies between qubit j and qubit j + k in the process of eliminating the non-nearest-neighbor interac- 
tions in preparation of cluster state (qubit j is any qubit in the molecule chain). The signs and numbers indicate the effective 
intense of the interaction between qubit j and qubit j + fc(Units: E+{a, kb)/4). 
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FIG. 1: Energy level structure of the quantum molecule system as a function of potential offset A between the three possible 
charge states. 
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FIG. 2: Relation between electron charge states and the interaction between two QD molecules, a is the distance between 
two QDs in one QD molecules and b is the distance between two QD molecules. Hollow circles stand for Quantum Dots, and 
solid circles stand for electrons. Every QD molecule(qubit) has two electrons. And the interaction of two qubit depends on the 
electrons distribution in the QD molecules. (a)(b) Effective interaction is zero, (c) Effective interaction is H+. (d) Effective 
interaction is H-. See text for details. 





FIG. 4: Drawing of two-dimensional QD molecule array. Each QD molecule is fabricated across the two layers of the system. 
The squares in the gate layer stand for the control gates, and the circles in the two-dimensional electron gas (2DEG) layers 
stand for quantum dots. When the gates are tuned properly, the corresponding qubit electron charge states can be modified. 



